Non-stationary parabolic equations for the 
quasi-monochromatic sound propagation in 
g ; media with a non-stationary background flow 



o 

(N 

Oh: 

00 



o 

O 



Ph. 



X 



M. Yu. Trofimov 

Il'ichev Pacific oceanological institute FEB RAS 

Baltiyskaya St. 43, Vladivostok, 690041, Russia 
e-mail: trofimov@poi.dvo.ru 



Abstract 



^ ^ The narrow and wide-angle parabolic equations for the quasi-monochromatic 

^»-^' sound wave packets propagating in a waveguide with a non-stationary 

, background flow are obtained. The results of numerical simulations 

are presented. 



in : 1 Introduction 

o 

Q ■ Parabolic equations for sound propagation in moving media were reported by 

. many authors, but any systematic derivation of such equations on the base 

O I of some universal asymptotic method was not presented in the literature 

I yet. In this paper we use the generalized multiple-scale method [T] to obtain 

some equations of this type. The material of this paper was presented at the 
seminar "Acoustics of inhomogeneous media X", Novosibirsk, Russia, 1-6 
June 2009. 



2 Derivation of equations 

We start with the equation [2] 

D — (—D p] - div — (- gradp^ +2 (-^ 
dt \pc'^ J dt \p J ^ dxj dxi \pdxj 

"^3 



where V = (Vi,V2,V3) = {u,w,v) is the velocity vector of the background 
flow written on the orthogonal coordinate system {x, y, z) (2;-axis is directed 
upward), Dt = {d/dt + ud/dx + wd/dy + vd/dz), p is the acoustic pressure. 

Let e be a small parameter, which is the ratio of characteristic wave length 
to the characteristic size of the medium inhomogeneties. Following the gener- 
alized multiple-scale metod [T] we rewrite ([T]) using the slow variables T = et, 
X = ex,Y = e^l^y, Z = t^l^z and the fast variable ^ = {l/e)e{X,Y, Z,T). 
We will assume that the flow velocity is 0(e) and then u = eU{X,Y, Z,T), 
w = e^/^W{X,Y,Z,T) and v = e^/^V{X,Y, Z,T), where the quantities U, 
W and V areO(l). We introduce for convenience the quantity n = 1/c and 
postulate the expansions 

= nl{X,Y,T) + euiX,Y,Z,T), 
p = po{X,Y,T,Z,r]) + ep^{X,Y,T,Z,r]) + ... , 

The density will be assumed depending on the slow variables only, p = 
p{X, Y, Z, T). At last, the partial derivatives are transformed by the rules 

dx dX dr] ' dy dY dr] ' 

analogously for the variables t z. In the rewritten in such a way ([1]) we collect 
terms with like degrees of e. The squares of the background flow velocities 
will be neglected as they was neglected in the derivation of ([T]). 

At 0(e^^ we obtain that the phase function 9 does not depend on variables 
Z Y,e = e{X,T). 

At 0(e° we obtain the Hamilton- Jacobi equation 

{exf - nl {Orf = {Ox + ^o^t) (^x - uoOt) = . (2) 

We shall consider the waves propagating in the positive X-direction and 
adopt as the Hamilton- Jacoby equation the first factor in ([2]). 

The equation at O(e^) does not contain pi on the strength of ([2]). Sub- 
stituting in this equation the anzats po = Ao{X,Y,T, Z) exp{iri), we obtain 
the equation for the amplitude Aq 

2inou- [Aq^x + uqAo^t] + ( -Ao,y ) + ( -^o,z ) + X^o = , (3) 
P \P Jy \P Jz 

where u = —9t is the local frequency of sound. 
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Assuming that 



where Cq = l/?^o; we obtain that 
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SO the potential of the parabohc equation contains the so called effective 
sound speed c + U, as was first shown in the paper [3]. We see that our 
equation is a generalization of the parabolic equation obtained in that work. 

Note that the dependence on variable Y in the obtained equation is just 
the same as on Z. Therefore in the sequel we shall not write the terms 
expressing the dependence on Y and consider the 2D waveguide. As will be 
easily seen, all result can be directly transferred to the 3D case. 

The analogous considerations at O(e^) gives the equation for the ampli- 
tude Ai, pi = Ai{X, Z, T) exp(i?7). 



2inoUJ- [Aix + uqAit] 
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xAi + J' = 0, 
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It can be shown [3] that the system of equations ([3]), ([5]) generalize the 
known wide-angle stationary parabolic equation obtained by the factoriza- 
tion method with the rational-linear Pade approximation of the square root 
operator. Moreover, even in the stationary case our equations contains the 
terms that the factorization method cannot produce. 



3 Initial-boundary value problems for the pa- 
rabolic equations 

For the simulation of the sound waves in the ocean the most interesting 
initial-boundary value problems for the Hamilton- Jacobi equatuion ([2]) and 
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the parabolic equations ([3]), © are the problems with X as the evolution 
variable in the domain Q = {Zq < Z < Zi} x {Tq < T < Ti}. 
For the energy norm 



E{X) -\9xAo\^ dZdT . 

Jto Jzo P 



under some simple and natural assumptions on the boundary conditions at 
dQ, the following theorem holds 

Theorem 1. The energy norm of the solution of the initial-boundary value 
problem for the parabolic equation ^ satisfies the inequality 

^(A-)<^(0).exp(£sup(^)..) (7) 

// uq can be represented in the form uq = no(X)a(T) then the following 
inequality holds 

E(X)<E(0)^supa(r). (8) 

This theorem immediately implies the uniqueness of solutions of such 
initial-boundary value problems for ([3]), (0) in functional spaces of the type 
Ci[0,X],L,m. 



4 Numerical simulation 

In the parabolic equations ([3]), dS]) the influence of the background flow is 
taken into account not only through the terms which explicitly contain the 
corresponding velocities, but also through the deformation of the sound speed 
and density profiles. As the geometric small parameter e in typical cases is 
much larger than the Mach number for the ocean medium, the question about 
essentiality of introducing the background flow velocities in the narrow-angle 
parabolic equation ([3]) is naturally risen. 

For an answer to this question the simulation of the harmonic sound wave 
propagation through the harmonic internal wave of a given mode was con- 
ducted, with the following parameters: the undisturbed waveguide of con- 
stant depth (100 m) was given by the parabolic sound speed profile with 
the minimum (1460 m/s) at the center of the waveguide and maximum 
(1500 m/s) at the top and bottom boundaries, the sound frequency was taken 
to be equal 100 Hz and the boundary conditions for the sound field was taken 
to be soft at the top boundary and hard at the bottom boundary. The point 
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Figure 1: Temporal variations of the transmission loss at the distance 6 km 
(left) and 12 km (right): with the velocity terms (dashed line), without these 
terms (solid line). The second mode of internal waves. 



sound source was situated at the center of the waveguide. The density strat- 
ification was given by the constant Brent- Vasala frequency (1/127 s~^), for 
the simulation were used the five minutes period internal waves of the second 
mode (wavelength=214 m) and the third mode (wavelength=142 m). 

The results are presented in fig. [T][^ in the form of the temporal varia- 
tions of the transmission loss at the distances of 6 km and 12 km from the 
source. We can conclude that the explicit introduction of the background 
flow velocities to the narrow-angle parabolic equation is essential. 



5 Conclusion 

In this paper the system of parabolic-like equations ([3]), (I5l), which can be 
used for numerical modelling of sound propagation in waveguides with non- 
stationary background flow, is obtained. Some properties of these equations 
are established. We hope that this information will be useful for the computa- 
tional acoustics community. Complete derivations, proofs and more extensive 
numerical modelling will be presented in the forthcoming paper. 
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Figure 2: Temporal variations of the transmission loss at the distance 6 km 
(left) and 12 km (right): with the velocity terms (dashed line), without these 
terms (sohd line). The third mode of internal waves. 
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